Causal Structure and Gravitational Waves in Brane World Cosmology 
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The causal structure of the flat brane universe of RSII type is re-investigated to clarify the bound- 
ary conditions for stochastic gravitational waves. In terms of the Gaussian normal coordinate of the 
brane, a singularity of the equation for gravitational waves appears in the bulk. We show that this 
singularity corresponds to the "seam singularity" which is a singular subspace on the brane universe. 
Based upon the causal structure, we discuss the boundary conditions for gravitational waves in the 
bulk. Introducing a null coordinate, we propose a numerical procedure to solve gravitational waves 
with appropriate boundary conditions and show some examples of our numerical results. This pro- 
cedure could be also applied in scalar type perturbations. The problem in the choice of the initial 
condition for gravitational waves is briefly discussed. 

PACS numbers: 04.50.-|-h, 98.80.Cq, 98.80.-k 



I. INTRODUCTION 

Since the proposal of a brane world model of our space- 
time by Randall and Sundrum [J (RS), the phenomenol- 
ogy of brane world cosmological models has been the sub- 
ject of intensive investigations in recent years. In these 
brane world models, our universe is regarded as a four 
dimensional boundary (brane) in a higher dimensional 
spacetime (bulk). Many authors have found more realis- 
tic models which include matter fields on the brane and 
realize the cosmic expansion 2], and tried to constrain 
on these models by the observational data. 

Due to the recent developments in the technology of 
astronomy, cosmological density perturbations and their 
observations through large scale structure and cosmic 
microwave anisotropics (CMB) have become the most 
stringent test to constrain on models beyond standard 
cosmologies. In order to examine the constraints on the 
brane cosmologies by the CMB observations, we have to 
know the 5-D informations about perturbations on the 
entire spacetime including the bulk 3ilJ|. 

In addition to the constraint by the CMB, the stochas- 
tic gravitational waves will be a promising candidate 
which provides a direct and even deeper test of such cos- 
mologies. One might see in principle earlier universe than 
the photon last scattering epochs by gravitational waves. 
Although equations for gravitational waves are simpler 
than those of density perturbations, the problem is es- 
sentially the same, i.e, one needs to clarify the evolution 
of gravitational waves not only on the brane but also in 
the bulk. Therefore, it is also instructive to clarify the 



minimum information to obtain the evolution of grav- 
itational waves before trying to clarify that of density 
perturbations in brane cosmologies. 

In order to give theoretical predictions of stochas- 
tic gravitational waves, many authors have adopted the 
Gaussian normal (GN) coordinate system in the neigh- 
borhood of the brane. In this GN coordinate system, the 
metric is given by the form 



(P{t, w) 



dr^ + (p{T,w)a^{T)dY.j^ +dw^, (1) 



and the equation of gravitational waves has the same 
form as that of the five-dimensional massless scalar field: 



□s/i = 0. 



(2) 
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The explicit forms of functions in the metric are given 
in the main text (see Sec. 1^). To obtain the theoreti- 
cal spectrum of stochastic gravitational waves, we just 
solve this equation with appropriate boundary condi- 
tions. However, Eq. ^ in terms of the coordinate system 
includes a singularity at w = Wh in the bulk, where 
the metric function (p vanishes. The treatment of Eq. (jJl 
near the singularity is one of difficulties when we obtain 
the evolution of stochastic gravitational waves Q- In 
fact, this singularity corresponds to the "seam singular- 
ity" discussed by Ishihara T\. 

The aim of this paper is to propose a numerical proce- 
dure to solve the evolution of cosmological gravitational 
waves avoiding the above difficulty in GN coordinate sys- 
tem. The procedure proposed here is based on the char- 
acteristic initial value problem according to the causal 
structure of the entire spacetime. This idea is analo- 
gous to the analysis of gravitational waves from a non- 
spherical domain wall by the one of the authors j^. We 
use a null coordinate instead of the proper time on the 
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brane. In this procedure, the boundary conditions in the 
bulk are replaced by the initial condition on a null hy- 
persurface and the above difficulty in the treatment of 
Eq. (121 near the singularity w = Wh is resolved if we sim- 
ply specify the initial spectrum on a null hypersurface. 

The organization of this paper is as follows: In Sec.lTTl 
we briefly review the global structure of a brane world 
universe and clarify the region covered by the GN coor- 
dinate in terms of the closed chart of the five-dimensional 
anti-de Sitter spacetime (AdSs). In Sec. IIIII we discuss 
the null hypersurface to clarify the causality of the prop- 
agation of gravitational waves in the bulk. In Sec. IIVI 
we develop the formulation to obtain the numerical solu- 
tions to Eq. (0) and show numerical examples of solutions 
which are derived by this formulation. The final section 
(Sec. 01 is devoted to summary and discussions. 

Throughout this paper, we consider the model with- 
out "dark radiation" following discussions in Ref.'^ and 
we only consider the flat Friedmann-Robertson- Walker 
(FRW) brane universe which is supported by recent pre- 
cise measurements of the CMB [lC|. 



II. GN COORDINATE SYSTEM IN AdSs BULK 

In this section, we first relate GN coordinate system to 
the flat chart in AdSs, and then consider the correspon- 
dence of GN coordinate system and the closed chart in 
AdSs which covers the entire AdSs. Through clarifying 
these relations, we can easily see that the region which 
is covered by GN coordinate system in the entire AdSs. 
Secondly, we find where the singularity of the equation 
(0 is in the bulk. 

Now, we consider a brane universe embedded in the 
AdSs with a negative cosmological constant As = —A/P. 
In terms of the static charts of AdSs , the metric on AdSs 
is given by 



K 



r\dY? 



(3) 



where K takes the values —1,0, and -t-1, corresponding 
to the negative, zero, and positive constant curvature of 
a maximally symmetric three dimensional space, respec- 
tively (for example, see Ref.Q)- The function jKiTK) in 
this metric is defined by 



When K = +1 and 0, the metric dT,^ is given by 



(4) 



dx+ + sin^ X+{dO'^ + sin^ 



'^^'^ " ^ dxl + xUdO' + sin^ 



{K = +l), 
(K - 0), 



respectively. Though we may consider the case K = —1 
which corresponds to the open FRW model, we do not 
treat this case in this paper. In these static charts, a 
trajectory of three brane is given by rx = rKitx) ■= 
a{T) and tfc = txiT), where r is the proper time of the 



world volume of the brane and a(r) is a cosmological scale 
factor on the brane. The equation of the brane motion 
is given by the generalized Friedmann equation 
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K 

7/2 



A4 

3 



^5 „2 



36' 



(5) 



where the dot denotes the derivative with respect to t, 
p is the energy density on the brane, A4 = nfX^ /12 + 
3As /4 is the cosmological constant induced on the brane, 
Gn = K|A/487r and ks a-re the four-dimensional and five- 
dimensional gravitational constants, respectively. Here 
and hereafter we assume Z2 symmetry across the brane. 

Cosmological solutions of the RS type brane world in 
terms of GN coordinate system were found by several 
authors 0], which is given by Eq. In this metric, 
functions ip and (p are given by 



(^(t, - J [(1 - ^)e^ + l)e- 



and 



iP{t, w) = (p{t, w) 



1 
2H 



dip 



(6) 



(7) 



where A := Vl + PH"^- We have chosen the coordinate 
w so that w = on the brane. As noted in Sec. HI this 
coordinate system has a coordinate singularity at w = 
Wh, where Wh is determined by the equation ip{T,Wh) = 
for each r. This equation yields 



exp I 



A + 1 
A-1' 



(8) 



This is nothing but the singularity in Eq. Q. 

To give the explicit coordinate transformations be- 
tween the metrics and (PJ with K — 0, we show the 
explicit forms of the functions (toi^o) in terms of (t, w). 
The function ro in Eq. ^ corresponds to the volume ele- 
ment on the T = const, hypersurface in the brane. Then, 
as shown in Ref. TT] , this function is given by 



(9) 



The cosmological scale factor a(T) is given by a(r) := 
rQ{T,w — 0). On the other hand, the explicit expression 
of to as a function of (t,w) is given by 



A + 1 



_ f 

A-l 



{I- A)e^f +1 + A 



, (10) 



where if,(T) is chosen so that to = tf,{T) on the brane 
{w = 0) for any r. The derivation of this expression H10(l 
is shown in Appendix ^ Clearly, the singularity in the 
equation |2Jl is just on the region ro = in the flat chart. 
Actually, in the vicinity of the singularity {w Wh), ro 
and to behave as 



ro 



a(r) A- I 



A+l 



1 



2(yl+l 



^0 — tb 



a{r) 



A + 1 



2 

A-IK-,^' 



o(a,(ii) 



(12) 
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where ^ := (1 — A)e'^T + 1 + ^4. To find where is the sin- 
gularity in Eq. ^ in the bulk, we first see the region of 
To = in the entire AdSs using the closed chart. On this 
closed chart, we can easily specify the point of the sin- 
gularity in Eq. (0) in AdSs bulk by tracing the spacelike 
geodesic which normal to the brane for each r. 

Now, we consider the relation between the flat chart 
and the closed chart of AdSs. Though this is already 
given by IshiharaQ, we repeat his arguments to find the 
event (p{T,Wh) = for each t in the closed chart. The 
AdSs is identified with the universal covering space of a 
hyperboloid, 




r*cosX+ 



FIG. 1: Embedded flat FRW brane. Hypersurfaces at which 
T or xo are constant are also shown in the flgure. 



(13) 



in six-dimensional flat spacetime with the metric, 

5 

ds^ = -dY^ - dY^ + dY^ . (14) 



In addition to these coordinate systems, it is conve- 
nient to introduce the temporal-radial subspace which 
makes the metric conformally flat 



/+(r+)(- 
arctan(r- 



-dtl 
70 



drl) 



t2 



(18) 
(19) 



The flat chart and the closed chart (Eq. cover this 
hyperboloid as follows: For the K = +1 case (the closed 
chart), 



/Yo\ 




Yi 




Y2 




Y3 




Yi 




\yJ 





+ rl cos{t+/l) 

T+ COS x+ 
r+ sin x+ cos 6 
r+ sin x+ sin 9 cos (p 
r+ sin x+ sin 9 sin 



(15) 



For the K = case (the flat chart). 



(Y2+Ya\ 




Y2-Ya 




Yi 




Yz 




Yi 




\ Y, ) 





roxo cos 6* 



V 



roXo smf/ cos < 
roXo sin 6* sine; 



(16) 



We should note that the flat chart does not cover the 
entire AdSs since ^2 + ^0 > (ro > 0). Comparing 
Eqs. I|15|) with Eqs. ptil) . we found the coordinate trans- 
formations from the flat chart to the closed ones: 



for the purpose of the investigation of the causal struc- 
ture. Note that in these coordinates, radial null rays are 
represented by straight lines at ±45 degree. In the co- 
ordinate system (i-|., cosx+, sinx+), the flat FRW 
brane is embedded as shown in Fig. ^ The r = const 
and Xo — const hypersurfaces on the flat FRW brane is 
also shown in the same figure, respectively. 

Through the relations ((Tn|l . and lfT7|l . we can see 
how r = const, hypersurfaces in the coordinate system 
foliate the entire AdSs bulk. To do this, we have 
to consider the spacelike geodesies normal to the brane 
specifying their starting point on the brane. Note that 
on each t = const, hypersurface in the brane {w — 0), 
both functions to(''") E^nd 7'o(r) are also constant, while 
the other coordinates on the brane are arbitrary. There- 
fore, we may choose a point (ioj^o) on the fiat chart as 
a starting point of these spacelike geodesies normal to 
the brane. Fixing the proper time t on the brane, the 
coordinate functions r^^ t^, and x+ behave as 



A 



air) _ 1 



I 

COSX+ 



arctan 
-1 



I ' 



(20) 



cosx^ 



smx4 



ra - 

2r+ 
roXo 
r+ 



1 



tan I y - 7re(io) 



t^ 



xl 



^2 



2 2 



I 

2to 



4 
P 



where 0(io) is the step function, i.e., 6(io) 
and e(to) = 1 for to > 0. 



(17) 



Xo 

for to < 



in the limit w Wh {ip{T,w) — > 0). This shows that the 
coordinate singularity at which ip[T,Wh) = is just on 
the "seam singularity" t+ —r^, cosx+ = — 1 (see Fig.|21l. 

Note that the seam singularity can be replaced by the 
other regular portion of a spacetime according to the cre- 
ation scenario of the brane universe. In this sense, we do 
not have to be afraid of this singularity, or equivalently, 
the singularity in Eq. seriously. Moreover, when we 
consider the characteristic initial value problem, the ini- 
tial conditions are chosen on a null hypersurface and the 
seam singularity is contained in its causal past. There- 
fore, the difficulties in the treatment of the singularity in 



4 



Eq. |5J are simply reduced to the choice of initial condi- 
tions for gravitational waves as seen in below. 



(t+/ 1) 



371 

T 




t+ = 

(seam singularity) 



r./l 



FIG. 2: Geodesies normal to the brane for each r (dot-dashed 
lines) in the entire AdSs spacetime. The geodesies converge 
to the seam singularity t+ — r* at which </p(r, w) = 0. 



III. NULL HYPERSURFACE AND CAUSALITY 
IN THE BULK 

To clarify the causal structure of the brane universe 
and the domain of dependence of gravitational waves, let 
us now consider the future directed radial null geodesies 
from the brane into the bulk. These null geodesies deter- 
mine the boundary associated with the causal region of 
the events on the brane. 

In terms of the flat chart, the equations of these null 
geodesies are given by 



dtp 



dX 



dxo 
dX 



0, 



(21) 



where E{> 0) is an integration constant and A is the 
affine parameter along the geodesies. Because we con- 
centrate on the radial null geodesies from the brane, xo 
is constant along the geodesies and the remaining equa- 
tions in Eqs. (|21|l yield 



Here, we note that the flat chart of AdSs has the fu- 
ture Cauchy horizon in the entire AdSs. This Cauchy 
horizon appears as the "infinity" {to,ro) (+oo,-|-0) in 
the flat chart. The above radial null geodesies from the 
brane approaches to this Cauchy horizon as the affine 
parameter A increases. We choose the affine parameter A 
so that A = corresponds to the point on the brane and 
A = > corresponds to that on the Cauchy horizon. 

On the other hand, in the closed chart, these null 
geodesies are given by 



2 1 

rl = — 
+ 4 



1 + - Xo Ko 



2C 



2 2 



C0SX+ = 



1 

2^4 



1 + ^ - Xo ) '^o + 2C 



(23) 



where C :— — ^ + ti,. In Eqs. 123|l . to and rg are param- 
eterized by A through Eqs. 1)22(1 . In the limit where the 
null geodesies from the brane approach to the Cauchy 
horizon (ioi^o) ^ (+oo,+0), Eqs. H23|l behave 



ICI, 



I 



arctan 



+ Tr, COSX+ = sgn(C). (24) 



These do represent a point on the future Cauchy horizon 
of the flat chart. 



(t+/i) 



371 

T 




, r)^=const 
(creatior^ 



t+= r* 




r*/l 



to 



PEX 



ni-EX + n) 



+ tb, ro 



-EX + n. 



(22) 



In Eqs. (|22ll . t\j and are constants which correspond 
to the initial condition for the null geodesies. We regard 
that the point (tb, rb) is on the brane in the flat chart and 
this point is the starting point of the geodesies. For this 
reason, we regard tb and rb as functions of the proper 
time T on the brane and these functions are constrained 
by the generalized Friedmann equation Eq. 



FIG. 3: Brane trajectory (lines) and future null geodesies from 
the brane (dot-dashed lines). Singularity where <^(r, w^) = 
is shown as a thick line, which is entirely contained in the 
shaded region 'D{tc). The evolution of gravitational waves is 
determined by boundary conditions on the brane and an ini- 
tial spectrum on dVijc). The initial spectrum is determined 
by the creation scenario of the brane universe. 

Actually, from the closed chart metric Q with K = 
we can easily see that the event represented by 



5 



Eqs. (|24|l is pointlike when we fix the proper time r on 
the brane. We also note the fact that Eqs. (|24|) do not 
depend on xo- This implies that all null geodesies with 
different xo from the same t — const, hypersurface on the 
brane reach to the same point represented by Eqs. (|24|1 . 
These null geodesies generate a null hypersurface. Hence 
we have one-to-one correspondence between a r = const. 
hypersurface on the brane and a null hypersurface in the 
bulk. We denote this null hypersurface associated with 
r = const, hypersurface on the brane by dT>{T). We also 
denote the causal past of d'D{T) by T>{t). We note that 
T>{t') C f (r) for any r' < r. In particular, the bulk 
region covered by the flat chart is foliated by the set of 
null hypersurfaces {9I?(t) | t is the proper time on the 
brane}. Since 9I?(t') C T>{t) for r' < t, the seam sin- 
gularity, which is also represented by dT>{T — s- +0), is 
entirely contained in I?(t) for any r. These situations 
are schematically depicted in Fig. O 

Inspecting this causal structure, we can formulate 
the characteristic initial value problem of gravitational 
waves. Suppose that the flat FRW brane universe is 
created at the instance r = Tc. As depicted in Fig. |31 
the initial state of gravitational waves in I'(tc) affects 
to their evolution only through the null hypersurface 
d'D{Tc). Once we specify the state of gravitational waves 
in the region 'D{tc) according to the creation scenario 
of the flat FRW brane, we can specify the spectrum of 
gravitational waves on d'D{Tc). Hence, to determine the 
evolution of the spectrum of gravitational waves after 
T = Tc, we only have to specify the spectrum of gravita- 
tional waves on dT>{Tc) as an initial condition and trace 
the evolution of them through Eq. |(5J as seen in the next 
section. 



troduce a null coordinate and rewrite down the equation 
of gravitational waves in a single null coordinate system 
(see Fig.m. 
Now, we introduce the function u by 

u = to 4 + - , (25) 

where (ij,,?*^) determine the zero point of this function 
u. The function u is constant on the null hypersurface 
d'D{T) for each t. We can easily confirm that the one- 
form {du)a is null through the metric © with K = 0. 




FIG. 4: Null coordinate and Gaussian normal coordinate in 
the fiat chart. 



IV. EVOLUTION OF GRAVITATIONAL WAVES 

At this point, we have a clear strategy to tackle the 
problem of the evolution of gravitational waves in the 
brane universe. The singularities ip{T,Wh) = in Eq. ^ 
are entirely contained in 2?(t), and the evolution of grav- 
itational waves are determined by an initial spectrum 
on its boundary 91? (r) and boundary conditions on the 
brane. Then, we do not have to care the singularity 
ip{T,Wh) =0 any more and there is no need to introduce 
any artificial boundaries in the bulk to impose boundary 
conditions at bulk infinity. In the following subsections, 
we first develope the formulation to obtain the numer- 
ical solutions to Eq. Q in the context of brane world 
cosmology (in Sec. IIV All , and then, 



we show numerical 
examples which are obtained by applying the formulation 
developed here (Sec. IIVB)l . 



Using the relations listed in the appendix fEa. (|A6l 
see also Ref. HtI), the derivative du is given by 



-(^] F'\T,w)dT-F-\T,w)dw 



where 



F{t,w) 



f[ro{T,w)] 



^f[ro{T,w)] +d'^ + a 



(27) 



In terms of the new coordinate system u and w ~ w, the 
five-dimensional metric on AdSs is written as 



ds' 



-F{T,w)du'^ - 2F{T,w)dudM + rl{T,w)dJ:l,i28) 



A. Numerical Formulation 

In order to develop the characteristic initial value prob- 
lem associated with the null hypersurface &D{t), we in- 



As mentioned in Sec. ^ the equation for gravitational 
waves in the bulk is simply given by that for the five- 
dimensional massless scalar field Eq. (jSJ. In terms of the 
coordinate system (m, w) Eq. |(5J) is given by 



6 



F 



-d—d 



3 



dro\ 
du , 



chu- 



drp 
dw 



dF 
dw 



3 



dro\ 
dw I , 



h{u,w;k) = , (29) 



where — fc^ is the eigen value of the Laplacian of dEo. 
Eq. (|29|l can be formally integrated so that 



duh{u, w) 



'(t,w) 

S{w')rl^^dw' + C'iu) ] , (30) 



where 



Siw) = 



F 
"2 



1 

2dW 



3F d In rg 
~2 9^ 



3 9 In To 
2 (9u 



^/i .(31) 



The quantities (9ro / 9u, Otq / dw and dw which appear 
in the definition H31|) are given by known functions as 



drp 
du 
drp 
dw 

dF 
dw 



dF 



drp 
dw 



drp 
dr 



d£ 
97 



dF\ 
dw I . 



(32) 
(33) 
,(34) 



These equations H3Uf) and H31|) are the main result of 
this paper. In these equations, initial data of gravita- 
tional waves are set by choosing the function of h{u,w) 
on dD{T) where the null coordinate u is constant on 
dVij). The function C{u) in Eq. (|30|) is determined by 
the boundary condition at the brane {w = 0) by 



duh{u,Q)^rp^''^{T,Q)C{u) 



(35) 



When the anisotropic stress on the brane is not induced 
due to matter fields on the brane, the boundary condition 
for gravitational waves h at the brane is the Neumann 
type j3| . This is accomplished by imposing 



duh{u, 0) = - ( ^ 1 

= F{T,0)d^h , 



(36) 



at the brane. Eq. H30() with the boundary condition l|36() 
are easier to solve than Eq. ^ and enough to predict the 
cosmological evolution of gravitational waves. Note that 
our method have an important advantage that applica- 
tion to high energy epochs is easy and straightforward. 
Actually, this advantage is seen in the numerical exam- 
ples shown in the following subsection (see Figs. [3 and 



B. Numerical Scheme and Examples 

Here, we show some numerical solutions to Eq. ^ 
in the context of the brane world cosmology. We first 
comment on some details of the numerical scheme to ob- 
tain solutions to Eq. H29|) with the boundary condition 
Eq. 

To solve Eq. H29|l . we evaluate Eq. H3U|I on each null 
hypersurface 92?(t(u, w = 0)), numerically. Once h{u,w) 
on dT>{T{u^w = 0)) is given as a function of we can 
evaluate h(u + du,w) on dD{T{u + du,w = 0)) by 

h{u + du, w) — h{u, w) + duduh{u, w). (37) 

Through this evaluation on each dT>{T{u, w = 0)),we ob- 
tain the gravitational waves h{u,w) in the causal future 
of the initial surface dT>{Ti„i) with an appropriate initial 
condition. By evaluating h(u, 0) on the brane, we see the 
behavior of gravitational waves on our brane universe. 

To obtain duh{u,w) by Eq. (|^ . we have to evaluate 
Eqs. on each dV{T{u,w = 0)) and Eq. ^ 

at the brane. To accomplish this, we evaluate a(T) — 
a{T)H{a{T)), F{t,w), and rp(T,w) on each d'D{T{u,w — 
0)). Since these functions depend on u only through r, 
we evaluate the change of the function t(u,w) on each 
d'D{T{u,w — 0)). Further, it is also convenient to use the 
scale factor a(r) at the brane as an time coordinate when 
we clarify the behavior of gravitational waves during the 
cosmic expansion. 

We evaluate the change of the function a{T(u,w)) on 
d'D{T{u,w — 0)) from Eq. by setting du = 0, which 
yields 



dw ) .„ " V^T j„ 



(38) 



where {drp/dT)^ is directly given by Eq. ©. By in- 
tegrating Eq.|(2Hl) with the boundary condition a(T) = 
rp{T,w — 0) at the brane, we obtain the function 
a{T{u,w)) along dT>{T{u,w = 0)). Since the relation be- 
tween a and r is given by the explicit integration of the 
generalized Friedmann equation Eq. (jSJ , we can evaluate 
F(t,w) and rp{T,w) on each dT>{T{u,w — 0)) by using 
Eqs. © and ijnji. The relation between the time step du 
and dr on the brane is given by 



dr — d 



^] F[r,w 



0)du , 



(39) 



which is led from Eq. H2t)|) by choosing dw — and w = 0. 
After proceeding the time step, we repeat the above eval- 
uation of air) and t{u,w) on d'D{T{u + du,w = 0). Thus, 
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we can evaluate F{t,w), ro(r, w), and d on each null hy- 
persurface d'D{T{u,w = 0)) from Eqs. ©, and (jSEJ, 
and hence, Eqs. (|32|l - 134|l . Through these evaluations, 
we obtain duh{u,w) by Eq. H3U|I and gravitational wave 
h{u, w) in the causal future of the initial surface d'D{Tini). 

When we carry out this numerical scheme, we have to 
specify the zero point of the function u. Suppose that 
we start our numerical calculation from the high energy 
epoch {PH^ ^ 1, or the dominated era). In this 
epoch, the generalized Friedmann equation and Eq. (|A3|) 
lead 

i/» ~ Ir-^ , (40) 
^^ [l4p-,.r)"' ,41) 

where p^Q and Oq are the energy density of radiations 
and the cosmic scale factor on the brane at the present 
universe, t'f^ j„j is the trivial initial value of ti, and we 
choose so that tj, jj^j = 0. Since the zero points of u and 
tb are arbitrary, we choose so that u = t'^ = and = 
iCTini) in Ea. (|25|l at the starting point of our numerical 
calculation r = Tini = l/4i?(aini). 

Following to these numerical scheme, we obtain numer- 
ical examples depicted in Figs. El and El In these exam- 
ples, we simply choose h = const, as initial conditions of 
h on the initial null surface d'D{Tini)- We take this initial 
condition only as an example to demonstrate the calcu- 
lation and investigate how the evolution of gravitational 
waves can be modified in high energy epochs qualita- 
tively. The wave number is chosen k — a(rini)_ff (Tjni), 
which corresponds to the mode just crossing the cosmic 
horizon at t = Tjni. In Figs. El and the scale factors 
of our examples are taken to be a(Tini) = 1.58 x 10"^^ 
and a(rini) = 1.58 x 10~^^, respectively. The expan- 
sion of the universe is completely dominated by radia- 
tion in such early epochs. We only consider ordinary 
radiation (CMB) and three species of massless neutri- 
nos for simplicity neglecting other matter components, 
and set the curvature scale of AdS5 bulk to be 1 mm. 
In this cosmology, the transition from the p^ dominated 
universe to the standard radiation dominated one occurs 
at a w 4 X lO^i"^. 

The evolution of gravitational waves in brane world 
differs that of standard cosmology mainly in two points. 
First, the amplitude of gravitational waves in brane world 
models becomes smaller than that in four dimensional 
standard models. This is due to the fact that decelerating 
brane motion excites the modes propagating into the bulk 
Second, the effective frequency lu of gravitational 
waves measured on the brane becomes larger due to their 
momentum not only along the three dimensional space (k 
in Ea. H29|) ') but also to the bulk direction. These effects 
becomes remarkable in the high energy epoch (PH^ ^ 



!)■ 
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FIG. 5: Evolution of gravitational waves on the brane in brane 
world cosmology (solid line) which cross the horizon in low en- 
ergy epochs [f'H^ <g 1). Also shown is standard 4D evolution 
of gravitational waves (dotted line). 
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FIG. 6: The same as Fig. El but for gravitational waves which 
cross the horizon in high energy epochs [l^H^ 3> 1). 



V. SUMMARY AND DISCUSSION 

In this paper, we have carefully investigated the causal 
structure of the flat FRW model of RS type II brane 
world and proposed the single null coordinate system to 
solve the cosmological evolution of gravitational waves. 
We have explicitly seen that in this null coordinate sys- 
tem, we do not have to care the singularity in Eq. 
Further, it is not necessary to introduce any artificial 
boundaries to impose some boundary conditions at the 
bulk infinity. Thus, we have shown that the problems of 
the singularity in Eq. ^ and the boundary conditions 
for gravitational waves at the bulk infinity are resolved 
if we simply choose an appropriate initial conditions for 
gravitational waves. 

The initial condition for gravitational waves which we 
considered to obtain the numerical examples (FigEl and 
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in the main text might not be reahstic one, since 
h = const, is not an exact solution to Ea. l|29|) but a 
mere approximated solution in the long wavelength limit 
k —^ 0. Since the aim of this paper is to propose a nu- 
merical procedure to solve the evolution of cosmological 
gravitational waves, the details and quantitative studies 
of numerical results and the problem on realistic initial 
conditions of stochastic gravitational waves are beyond 
the current scope of this paper and they will be investi- 
gated in our forthcoming paper However, it is in- 
teresting to discuss the evolution of gravitational waves 
with an appropriate initial conditions and clarify the fi- 
nal spectrum of stochastic gravitational waves resulting 
from various creation scenarios of the brane world. 

The initial conditions for gravitational waves in brane 
world cosmologies crucially depend on the creation sce- 
nario of the FRW brane. Many kinds of cosmological 
scenarios have been proposed so far ^T^. Among them, 
there are some scenarios where the FRW brane is created 
after the inflationary phase. If we adopt these brane in- 
flationary scenarios, it might be natural to consider that 
the initial spectrum of gravitational waves is determined 
in this inflationary phase. The spectrum of gravitational 
waves under the deSitter evolution of the brane is dis- 
cussed by several authors 0, 0, 0. It was pointed 
out that gravitational waves decay away except but the 
zero-mode, and it approaches to a constant amplitude in 
the inflationary phase Q. These discussions are based 
on the GN coordinate system. Since GN coordinate sys- 
tem does not cover the entire bulk space, these discus- 
sions seem inappropriate to specify the initial spectrum 
of gravitational waves. However, it was also shown that 
the vacuum defined on the deSitter slicing asymptoti- 
cally approaches to the vacuum defined in terms of the 
Poincare coordinates on AdS^ • This will imply that 
the vacuum state on the static AdSs frame is appropriate 
as the bulk initial spectrum of gravitational waves when 
we consider these inflationary scenarios. 

If we do not adopt the inflationary scenarios, we have 
to specify the initial spectrum according to the other cre- 
ation scenario of the FRW brane. However, in any case, 
once given the initial conflguration on a null hypersur- 
face, our method can be applied to solve the evolution of 
gravitational waves. The final spectrum of the stochastic 
gravitational waves can be a powerful probe to inves- 
tigate the existence of extra-dimensions by comparison 
with the spectrum in the four-dimensional standard cos- 
mology. We leave this to future works. 

Besides the evolution of stochastic gravitational waves, 
the procedure developed here will be applicable to discuss 
the evolution of the density perturbations in the brane 
world. The problem in the choice of the initial spectrum 
of the density perturbation will also arise as discussed 
above and this initial spectrum will also depend on the 
creation scenario of the FRW brane. However, accord- 
ing to the causality discussed in this paper, we can easily 
expect that the density perturbation is completely de- 
termined by the initial condition on a null hypersurface 



and the boundary conditions at the brane. Though this 
expectation should be confirmed by examining the equa- 
tions for the density perturbations of the brane world, it 
is quite interesting to compare the evolution of the den- 
sity perturbations in the brane world scenario with that 
in the conventional four-dimensional cosmology. We also 
leave this problem as a future work. 
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APPENDIX A: GAUSSIAN NORMAL 
GEODESICS 



We briefly review the relation between the Gaussian 
normal (GN) coordinates and the flat chart in AdSs. To 
do this, we flrst consider the spacelike geodesies normal 
to the brane. Let us consider the unit normal vector 
«° {gabn°"n!' = 1) tangent to these spacelike geodesies. 
The existence of the Killing field {d/dto)'^ in the bulk 
spacetime ensures that the existence of the constant £ = 
— 5afcf^"(9/9^o)'' along the geodesies. Using this constant, 
the components of the normal vector n° are given by 



d_ 
£ f d 



d 
dro 



firo) \dto 



VfW) 



d 

dro 



.(Al) 



The explicit orbit of these geodesies is given by the 
integration of dx^ = n'^dw. For tq component, we obtain 



2r' 



VE^ = VE^ cosh 



--[w ~ Wo) 



(A2) 



where wq is the additional constant of integration. We 
choose E and wo by the conditions: (i) the geodesies are 
normal to the brane world volume of tq = R{to) = air) 
at to = ti, (i.e. Ua oc [ro — R{to)])', (ii) w ~ on the 
hypersurface. Since dto/dr along the brane is given by 
the definition of the cosmic time t: 



dr' 



-f{ro)dtl 



-\ro)drl 
'\ro)d'dT' 



(A3) 
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the above two conditions lead 



Wo{tb 



I 



■ cosh 



2a^ + P£^{tb 



(A4) 



Thus, substituting eas. (|A4|) into ea. (|A2|) . we obtain eq. 
. On the other hand, the integration of the tg compo- 
nent of dx^ = n^dw, i.e.. 



dtp 
dw 



n*0 



Pair) 



(P{t, w)a^(T) 



(A5) 



leads Eq. |Tn|). 

Finally, we summarize some useful relations. The fact 
that T is constant along these spacelike geodesies leads 



dto 
dw 



firo) ' 



(A6) 
(A7) 



Also, from the definition that dw = gf^i,n^dx'^ , we obtain 



dw 
dw \ 



a(r) , 

Vfiro)+d^T) 
firo) 



(A8) 
(A9) 



These relations are used in the main text (see Sec lIV|) . 
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